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Abstract 

We determine how the yield of large transverse momentum hadrons is modified due to 
induced gluon radiation off a hard parton traversing a QCD medium. The quenching factor 
is formally a collinear- and infrared-safe quantity and can be treated perturbatively. In spite 
of that, in the pj_ region of practical interest, its value turns out to be extremely sensitive to 
large distances and can be used to unravel the properties of dense quark-gluon final states 
produced in heavy ion collisions. We also find that the standard modelling of quenching by 
shifting p± in the hard parton cross section by the mean energy loss is inadequate. 



1 Introduction 



The so-called jet quenching 1T|-§| is considered an important signal of the production of a new 
state of dense matter (quark-gluon plasma) in ultrarelativistic heavy ion collisions. This is 
understood as the suppression of the yield of large transverse momentum jets or particles with 
respect to proton-proton collisions. 

In this paper we concentrate on the quenching effect in inclusive particle spectra, due to the 
energy loss by medium induced gluon radiation 0-^. 
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Inclusive production of particles with large p±, in, say, proton-proton collisions can be 
parametrized as a power 

' n = n{p±) = -— In (1) 

dp± P± dmp± dp]_ 

with n an effective exponent which slowly decreases with increase of p^. In reality, for moderately 
large p± the exponent in (P is much larger than the asymptotic value n = A corresponding to 
the limit p^ ^ oo {s/p]_ = const). 

The effects that add to the value of n are: 



1. x-dependence of the parton distributions which decrease with increase of parton energies, 

Xi,X2 OC \/p\/s, 

2. the bias effect due to vetoing accompanying gluon radiation off the primary large-pt parton 
that produces the triggered particle with p_L < pt, 

3. running of the coupling the parton-parton scattering cross section. 

As a result, in practice n is seen to be as large as 10. In what follows we shall treat n as a large 
numerical parameter neglecting relative corrections of the order of (Such an approximation, 
though unnecessary, allows us to derive a simple analytic expression for jet quenching due to 
medium effects.) 

In the presence of a medium the inclusive spectrum changes. Multiple interactions of 
partons with the medium lead to two competing effects. One the one hand, multiple scattering 
in the initial (as well as in the final) state partially transforms longitudinal parton motion into 
transverse one, thus enhancing the yield of large-p^ particles (the so-called Cronin effect [p5| ). 
On the other hand, medium induced gluon radiation accompanying multiple scattering causes 
parton energy loss and therefore suppresses the particle yield. 

At large transverse momenta the second effect takes over. 

To find the inclusive particle spectrum, one has to convolute the production cross section of 
the parton with energy p± + e with the distribution D{e) in the parton energy loss e in the final 
state:Q 

^^mcdium(p^) _ f ^ , da~^\p^ + e) 

dpi -j^'^^'^ dpi ■ ^' 

Since the effects of the standard jet fragmentation which reduce the energy of the leading particle 
in the jet are already included in the vacuum cross section, the distribution -D(e) here describes 
specifically the additional energy loss due to medium induced gluon radiation in the final state. 

*For the sake of simplicity we consider particle production at 90° and equate the transverse momentum of the 
particle p± with its energy. 
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The quenching effect is customarily modelled by the substitution 



(3) 



The shift parameter S in is usually taken either proportional to the size of the medium, 

S = const ■ L , (4a) 
or equal to the mean medium induced energy loss [|ri[] 

S = AE= [ dttD{e) oc a, . (4b) 



The former ansatz has no theoretical justification while the latter emerges as a result of the 
Taylor expansion of (0) based on the e <^ approximation: 

deDie) . ^^iPl±A = f deD(e) ■ ^+ / deeD(e^ ^ ^^^^^^ ■ 



dpj_ J dpj_ J dp± dpj_ 

(5) 

da d f da \ da{p± + AE) 

' AE ■ 



dp]_ dpj_ \dpj_J dpj_ 

Such an approximation misses, however, one essential point, namely that the vacuum distribution 
is a sharply falling function of p±. This causes a strong bias which leads to an additional 
suppression of real gluon radiation. As a result, the typical energy carried by accompanying 
gluons turns out to be much smaller than the mean ( |¥5D . 

In this paper we study an interplay between the energy loss and the cross section fall-off and 
show that it leads, in the region of transverse momenta of practical interest, to the pj_ dependent 
expression for the shift 

S{p±) oc ^/pZ- 



2 Medium induced energy loss 
2.1 Transport coefficient 

The two effects — parton transverse momentum broadening and medium induced radiation are 
closely related and are determined by the so-called "transport coefficient" q which characterizes 
the "scattering power" of the medium [^: 

g(^) = p/cig^g^^. (6) 

Here p is the density of scattering centres, and da^^^ is the single scattering cross section for 
a projectile parton in the colour representation R, with Cr the corresponding colour factor 
(Cf = {Nl - l)/2N^ = 4/3, Ca = Nc = 3 for quark and gluon, respectively). 
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The dimensionless ratio q/ p characterizes the "opacity" of the medium for an energetic gluon. 
It includes the region of small momentum transfers where the perturbative treatment is hardly 
applicable, and should be regarded as (the only) unknown medium-dependent parameter of the 
problem. 

In "cold" nuclear matter q can be calculated perturbatively and related to the gluon density 
[xG^x^Q"^)] of the nucleus at a low momentum scale Q"^ '::=l qL and small but not too small x, 
where the gluon density has little dependence on x [|lT|: 



(7) 

Hereafter we shall label q the gluon transport coefficient. Taking p = 0.16 fm~^, = 0.5 and 
xG{x) = 1 in (|^) results in 

~ 0.009 GeV^ ~ 0.045^^. (8) 

fm 

q enters, in particular, as the proportionality factor between an accumulated parton transverse 
momentum squared and the size of the medium traversed, oc qL. The quark transport 
coefficients extracted from experimentally measured transverse momentum nuclear broadening 
of Drell-Yan lepton pairs |]26| agrees with the theoretical estimate (IS). 



In the case of heavy ion collisions, the scattered hard parton traverses a medium that is 
expected to have an energy density much higher than that of nuclear matter, and the corre- 
sponding transport coefficient ghot can be much larger. If hot matter is formed in the final state, 
a perturbative estimate for the QGP with T = 250 MeV gives |]TT[] 



ghot ^ 0.2 GeV^ ^ 1^. (9) 

fm 



2.2 Induced gluon radiation 

To find the distribution D{e) in the parton energy loss e in the final state we need to recall the 
basic properties of gluon radiation caused by multiple parton scattering in the medium. 



Introducing the characteristic gluon frequency 



UJr 



(10) 



with L the length of the medium, the inclusive energy spectrum of medium induced soft gluon 



radiation {uj <^ p±) reads 11121,113 



a 

OJ 



In 



cos 



1 + 



2uj 



a 



In 



cosh 



2uj 



sm 



2uj 



2a,C 



a = 



71 
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This distribution peaks at small gluon energies, 



2^c 



, a; < cuc , 



(12) 



while for energies above the characteristic scale it is small and falling fast with uj: 

dl{ijj) a /cUc^ 2 



u- 



UJ > UJr 



(13) 



duo 12 \ 

The multiplicity of gluons with energies larger than a given uj is given by the integal of the 



inclusive gluon spectrum (|TT 



N{uj)= I duo' ) / = a I — In (cosh^ 2 - sm^ 2) . 



duo' 







z 



Evaluating this integral for x = cj/tUc <^ 1 we have (cf. (IT^)) 



N{uo) ~ a 



+ In 2 In X- 1.44136 + O (e-x.^{-^/2j^) 



with the constant term found by numerical integration of the exact spectrum ([TT|) . 



(14) 



(15) 



2.3 Energy loss distribution 

Since the vacuum spectrum (|l]) is falling fast with p^, the bias effect forces the distribution -D(e) 
to the smallest energy losses possible, e/p_L ^ 1. In these circumstances the "final" parton is the 
one that had been produced in the hard interaction: the quark-gluon transition is additionally 
suppressed as e/p_[_ and can be neglected. Bearing this in mind, we will treat -D(e) as the 
distribution in energy that the hard parton (a quark or a gluon) loses to medium induced gluon 
bremsstrahlung. 

The spectrum of the leading particle can be characterised by the probability -D(e) that the 
radiated gluons carry altogether a given energy e. The corresponding expression based on inde- 
pendent emission of soft primary gluons reads 



Die) 



00 ^ 

y- 

n=0 



W / duo 



dl{uOi 
duo 



uji ■ exp 



duo 



dl 

duo 



(16) 



where the last factor accounts for virtual effects. In a standard way, the energy constraint can 
be factorized using the Mellin representation, 



1=1 



c 



— e--n- 



(17) 



1=1 



5 



after which multiple gluon radiation exponentiates, and the answer can be written as 



D{e) 
DM 



'c 
exp 



duo ^ 



(18a) 
(18b) 



The contour C in (|T^ and ( |18a| ) runs parallel to the imaginary axis in the complex i^-plane, 
Rez/ = const. An equivalent expression can be written in terms of the integrated gluon multi- 
plicity ([T^) . Integrating ( 18b ) by parts, we obtain an elegant formula 



Diu) 



exp 
exp 



dz e~'N\- 



(19a) 
(19b) 



As we shall see shortly, the value of the dimensional variable vio^ is linked with the characteristic 
parameter of the problem nujc/pi_ which is typically much larger than unity. Therefore, in the 
essential integration region in ( p,8b| ) we have uj ~ 1/v ~ P:L/n <^ Uc, and the approximate 
expression (plSf ) can be used. This gives 



D{u) ~ exp {-a (V27rz/cjc - In 2 ln(z/tUc) - 1.84146)} 



(20) 



which expression allows one to evaluate the energy loss distribution D{e) for e ^ cj^ analytically 
in terms of the hypergeometric function. 

For the purpose of illustration we present here only a rough estimate based on the leading 
small-energy behaviour of the gluon distribution udl/doo oc ag^/uj^Juj which translates into 



Z)(z/) ~ exp |— a^/27u^} . 
The Laplace integral ( |18a| ) then becomes Gaussian and yields 



eD(e) 



a\l — exp 



-KO^UJc 

2^ 



(21) 



(22) 



This approximation is fine for illustrative purposes, in particular because the distribution (|2^) 
is properly normalized to unity due to -D(i^ = 0) = 1 in (^Tj) . As expected, in the first order 
in as the energy loss spectrum coincides with the probability of emission of one gluon with the 
energy = e, see (P^). The one-gluon approximation spectacularly fails, however, at small e. If 
the energy loss is taken as small as e ~ a^cUc, the distribution reaches its maximum at 

e = Tia^ ujc <^ oJc (23) 

and becomes exponentially small at smaller energies, due to form factor suppression. 
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3 Quenching 



We introduce the medium dependent quenching factor Q to represent the inclusive particle 
spectrum as 

= ^'^^ 

with 



<3(Px)= (25) 



Using the fact that the effective exponent n in the ratio 



R 



entering the convolution (^) is numerically large, we can replace R by the exponential form 



T16 

R = exp ( ) ■ 

P±J 



l + O 



2n 



(26) 



The accuracy of such a substitution is rather good all over the region of practical interest where 
the quenching suppression is moderately large. 



The exponential approximation (|2g) results in a simple expression for the quenching, as the 
suppression factor equals a given Mellin moment in e of the spectrum D{e) (see ( |18a| ) and (|I^)): 

Q(p^) ~^"rfeD(e) exp|--^-e| = ^ (^^) = {~ ^ (^^) } • (^7) 

Representing the quenching factor as 

Qip±) = exp -5(^^)1, (28a) 

S(p±) = -^\nD(—]= r duNiuj) exp (- — 1 , (28b) 
n \P±J Jo [ P± ) 

and using again the large-n approximation, we can cast the quenching effect as a shift of the 
vacuum spectrum: 

dpi - dpi ■ ^ ^ 

It should be noted that the shift approximation (|29|) underestimates suppression if the quenching 
factor Q becomes very small. 
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The general expression (^) actually solves the quenching problem, given the gluon emission 
spectrum dl/du. To verify the validity of the approach and the approximations made, we need, 
however, to estimate the characteristic gluon energies that determine the answer. To this end 
we proceed with a qualitative analysis of QCD quenching. 

Using the small-energy approximation ([T5|) in (^) we derive 

Q{p^) = exp{-v^iv(^)} ^ exp{-iv(|^)}. (30) 

This expression has a simple physical interpretation. It is a typical exponential form factor 
suppression whose exponent equals the probability of gluon radiation in the forbidden kinematical 
region or, in other words, the mean multiplicity of virtual gluons. We see that the energies of 
real gluons are cut from above as 

00 < , OOi = , [6l) 

IT n 

which (modulo a constant^ factor vr) equals a small (1/n) fraction of the energy of the registered 
particle, due to the bias effect. The characteristic energy parameter of the quenching problem 
can be directly expressed in terms of the hard vacuum cross section by the general relation 



^i(Pj 



d , da''™'' ' 
m 



dp_[_ dp]_ 

which does not rely on the power approximation (|T]) with n =const. 
For the "shift" function ( P^ in the approximation (|T3p we obtain 



(32) 



S(Px) - J!^^!^. (33) 
V n 

This result namely, the p±_ dependent shift, is very different from the models @) that were 
discussed in the literature and are being used to predict/describe quenching To understand 

the origin of (|33D in the rest of this section we consider and compare the mean energy loss with 
a typical energy loss that characterizes quenching. 



Mean energy loss. The mean medium induced energy loss of a parton with a given energy 
E is determined by the integral 



°° d ~ 

deeD{e) = - ^D{v) 



u=0 



duj ijj 



dl{uj) 
duo 



duN{u}) 



(34) 



where we have used (0). Given oodl/duo oc A^(ci;) oc uo this integral is determined by the 

(35) 



largest available gluon energies, oo < cUmax, resulting in [jlT 

^max 



Ai? oc as\JijOc UOn 



min { cUc , E} . 



■'dimensional, some would add, 



8 



The commonly used identification of tlie sliift with the mean energy loss (IB) is valid in the 
region of large particle energies p±_ > nujc, corresponding to the quenching parameter ^ 1. 
Indeed, since the integrated gluon multiplicity N{lj) vanishes fast for u > Uc, for z/ < l/tUc we 
can omit the exponential factor in the integrand of ( |19a| ) to derive 



S{P±) 



duoN{uo) = AE. 



(36) 



In this region, however, the quenching itself is vanishingly weak since AE oc agUJc and 



- In Q{p J 



AE oc a. 



P± 



< a.. 



Gluons responsible for the mean energy loss are rare, O (a^). Strictly speaking, rare fluctuations 
with energetic gluons (u; > ujc) in dl/duj do contribute to quenching via the virtual suppression 
factor in ( p.8a|) : 

dl{uj) 



SlnQ 



dto 



duj 



This contribution, however, is bound to be small (negligible) as long as the mean multiplicity 
(not mean energy^}) of such gluons is of the order as- 



This argument applies both to the "canonical" AE and to the additional contribution 

to AE originating from emission of energetic gluons when there is only a single scattering in the 
medium, found in a series of recent papers initiated by Gyulassy, Levai and Vitev pO|-p^. The 
GLV energy spectrum has an enhanced high-energy tail. 



; OC a , UJ > UJc 

duj ijj 



cf. 



^j{BDMPS)( 



duj 



oc a ^ — 



thus inducing a potentially large contribution to AE from the region uj ^ uj^- Such fluctuations, 
however, do not affect quenching since 



duj 



duj 



Typical energy loss. To have a significant quenching, — ln(5(p±) = 0{1), we have to have 
N{p±/n) ~ 1, according to (^71) , which translates into p± < a'^nuc <^ nuc- 

The value of the convolution integral @ that determines quenching results in an interplay of 
the steep fall-off of the parton cross section with e and the form factor suppression of small losses. 
To estimate the characteristic energy loss for a given p± we invoke the approximate expression 
for the distribution D(e) to write 



^, , a f°° dx f TTtt^ nuJc 1 

Q{p±) ^ ^ — exp<^-- x^=exp 

V2Jo X2 I 2a; p± } 




2p_L u, 



X = — . (37a) 
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In the region where the quenching is strong, Q <^ 1, the steepest descent evaluation apphes, 
provided the exponent on the r.h.s. of ( |37a] ) is large. 

With p± increasing, quenching becomes week, 1 — Q{p±) <^ 1. In this kinematical region we 
use the fact that D{e) is normalized to unity to write an equivalent representation 



1-Q{PA_) ^ \ — / ^■exp<^ — \, y = —. 37b 

\l ^p± Jo y y/y I ^p±y J p± 

Now, if the exponent in the last factor is small it can be dropped, and the integral over y is 
determined by y = O (1). 

Combining the two estimates we arrive at 



V — 2n — p± < -a nujc, (38a) 

(e) ^ — for p±>—a'^nuJc, (38b) 



where (e) is the typical energy value dominating the integrals (^7^. The two expressions (|3l 
match at the border value. 



Typical energy loss and the shift. In the first regime ( P8a| ) gluon radiation is omnipresent. 
Invoking (^) we observe that the characteristic energy loss (e) in ( |38a| ) equals half of the shift 
function S{p±). The factor | may look anti-intuitive at a first glance. To appreciate it we need 
to recall that the substitution 



da 



medium 



iP±) 



^(P± + (e)) 



dp\ 



dp] 



accounts only for a part of quenching: the second ingredient of the convolution (0) — the energy 
loss distribution D {{e)) — also supplies an (equal) exponential suppression factor, giving 



S{PA 



2(6) 



na^ujc p± 



2n 



(39) 



In the complementary region of larger pj_, the situation changes. Here gluon radiation is rare 
and can be treated as a correction. In the one-gluon approximation the loss is equal the gluon 
energy, (e) = Ui in ( |38bD . It has to be multiplied by the (small) gluon probability (multiplicity) 



AM(cu) 



const 



LU- 



dl{uj) 
dui 



(40a) 



to obtain the shift (B3| 



S{PA 



uoi ■ /^M{uJi) 



(40b) 
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Characteristic single gluon energy. In (plD we have introduced the characteristic (single 
gluon) energy loi as a border separating real and virtual emissions. Real and virtual gluons 
softer than Ui cancelled, while the multiplicity of harder virtual gluons provided the form factor 
suppression which resulted in quenching according to (0). Given an interpretation of uji as the 
maximal energy of real gluons, the shift can again be expressed as in ( [40b|) which evaluation is 
now valid for arbitrary p±, across (^). 



4 Approximations 

We start the discussion of the approximations made by noting that having written the factorised 
convolution we tacitly supposed that the leading quark turns into the registered hadron 
(hadronizes) outside the medium, which implies 

Then, the power approximation for the vacuum spectrum with n =const can be justified. Our 
derivation could have been damaged if the effective exponent n in (|l|) changed significantly over 
the energy range p± pj_ + e in the essential region in e. Substituting the typical energy loss 
( p8a|) we observe that it provides a correction 

which is actually relatively small for all values of Q of practical interest. This estimate allows 
us to evaluate n at the value of the transverse momentum of the registered particle, n = n{p±), 
and use the general expression (|3^ ) as the definition of the characteristic energy parameter Ui. 



Soft approximation. The medium induced gluon distribution (|Tl]) is valid in the soft gluon 
approximation, uj <^ p±. This approximation is justified by the fact that the main contribution 
to ( PUD comes from gluons with energies 

P-L 

u ~ Ui = — 'ti p± . (41) 
n 

Moreover, finite quenching, \nQ~^ = O (1), implies p± < nuc so that 

u ~ uji = — < Uc <^ ujc , (42) 
n 

so that the approximate expressions (|12|), (p!5[) are legitimate.^ 

^Subleading corrections due to hard gluons will be discussed in the next section. 
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Independent gluon radiation. The master equation (p!6[) was based on the independent 
gluon radiation picture. To estimate significance of possible interference effects we need to look 
at a typical density of gluons in configuration space. 

To this end we recall that the lifetime (formation time) t of a medium induced gluon with a 
given energy uj and transverse momentum fc^ follows from the relations 

t-^, klr^\^^\ (43) 

where is a typical transverse momentum transfer in a single scattering and A the gluon mean 
free path. This gives (/i^/A ~ q) 

(44a) 

Q 

kj_ ~ a/cJ^. (44b) 



Multiplying ( |44a| ) by the multiplicity density ( [40aj ) and dividing by the size of the medium, we 
estimate the "gluon occupation number" , 



Occup ~ t ■ AM ■ — ~ . — ■ a\ ■ — = a . 

L \ q \ u L 

This looks natural: a bare quark produced in a hard interaction builds up its wave function by 
emitting gluons, with a being the gluon density per unit phase space. The role of a (dense) 
medium is to strip the quark off these gluons, so that the formation process repeats again and 
again, resulting in induced production of many gluons. 

The smallness of the overlap between gluonsQ allows us to neglect possible interference effects 
in multiple gluon radiation as potentially contributing at the level of 0{as), while the main 
effect of the resummed medium induced radiation in (|16D is O (1). This justifies, aposteriori, the 
Poisson approximation (|T6]). 



Scale of the coupling. According to ( [44b| ) the characteristic scale of the QCD coupling 
as{k±) in the gluon emission spectrum (pHl) is rather small. For typical gluon energies (0) we 
have (g = 0.2 GeV^) 

*,x~(|p.)"* - 700 MeV.(^)"' =0(1 GeV), 

which value is rather low and marginally increases with p± reaching ki± ~ 1.5 GeV for p± ~ 
100 GeV (n ^ 6). In what follows we take a fixed value for the coupling 

a, = 0.5 , a = — ~ 0.42 , 

71 

motivated by the studies of the effective QCD interaction strength in the infrared region, 
similar argument applies to an interference between medium- and uacuMm-produccd gluons. 
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Bethe— Heitler limit. The perturbative approach to the problem inevitably limits gluon en- 



ergies from below. In particular, the gluon spectrum ( [Tl| ) describes the Landau-Pomeranchuk- 
Migdal (LPM) suppression of the independent Bethe-Heitler radiation and only applies to gluons 
with lifetimes larger than the mean free path, t > X, which is equivalent to limiting gluon energies 
from above, see fHF 



uj > cubh ^ jJ' X ~ . (45) 

In spite of the fact that the perturbative answer for the quenching factor (^) is formally infrared 
safe, its sensitivity to the region of small gluon energies may be significant, since the characteristic 
energy scale of the problem is strongly reduced by the bias effect. Indeed, according to (^), to 
keep the value of the quenching factor fully under perturbative control we have to choose the 
transverse momentum well above 

p± > nnuBH ~ 10 GeV, (46) 

where for the sake of estimate we have taken the vacuum exponent [0, n ~ 12 ^ 13 and 
cubh — 300 MeV corresponding to a hot medium with q ~ 0.2 GeV^ and A ~ | fm. At gluon 
energies comparable to cubh Debye screening effect should also become important thus limiting 
our ability of doing reliable calculations. 

To estimate the absolute minimal value of the quenching factor we take p± < uubh and, 
keeping in mind that N{uj) flattens out at small energies, from ( pT)) obtain 



In — - — ~ N{lubr) ' 



^ , ^ rs-/ 

[Q{P±)]mm ^'^ '^"^ TT V ^^BH TT A^ A/j ' 

where we have used 2uc/ujbh — {L/Xy and introduced the mean free path for the parton R 
(quark in our case) which is related with the gluon mean free path by CrXr = NcX. 



5 Illustrations 



5.1 Shift function and subleading (hard gluon) effects 

If the integrated gluon multiplicity in (^) depended on a single variable u/lJc, the quenching 
factor would have been a function of a dimensionless ratio X = In reality, such a scaling holds 
only approximately. Indeed, the gluon radiation spectrum ( [TTl) and, therefore, the multiplicity 
A^, depend on two energy ratios, lj/lUc and x = lj/p±. The latter takes care, in particular, of the 
phase space restriction^ on the energy of radiated gluons, dl / du oc N oc Q{1 — x) . 

According to (0), the main contribution to quenching originates from the region x ~ uJi/p± ~ 
1/n <^ 1. The correction coming from "hard" gluons with x < 1 turns out, however, to be 

^here we use an approximate equality of p± to tlie total jet energy 
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relatively large, O (l/y/n). To take into account non-soft corrections, in the numerical evaluation 
we supply the spectrum (|11|) with a factor (1 — x)Q{l — x) as given in [jl2|,[T^ 



The effect of the hard gluon correction is illustrated in Fig. ^ where we plot the shift normal- 
ized, for convenience, by AE = jaUc (the BDMPS mean energy loss |11]). 



Normalized shift 



S{P±) 
AE 



. 8 



. 6 



0.4- 



0.2 1 




0.2 



0.4 



. 6 



0.8 



1 P± 



Figure 1: Shift S{p±) as a function of X = ^ for n = 4 (lower curve) and n = 10 (upper curve). 

These subleading corrections can be treated analytically. With account of the phase space 
restriction, u < p±, and the "hard" factor (1 — x) in ([TT|) the integrated multiplicity (|1^) gets 
modified as 



X 



ln2 [Inx + 1- x] \ , 



X 



P± 



(47) 



{N{uj > p^) = 0) . 

The energy integral for the shift function ( P8b| ) then becomes 



Sip±) = p± dx N {xp±) e" 

'0 



2tt 



P± 
nojr 



1 



'Wn 2n 



P± 



In 2 



Inn + 7s - 1 + - 

n 



(48) 



where we have omitted exponentially small terms O (exp(— n)). We see that the hard correction 
is rather large, O {l/^/n), and significantly modifies the behaviour of the shift even in the strong 
quenching limit, X = — ^ 1. 



The exact numerical evaluation of S{p±) is compared with the approximate formula ( ]4^ ) in 
Fig. ^ for n = 4. 

The comparison shows that even for the smallest n value the approximation based on <^ Uc 
turns out to be rather good up to p_L = 2 ^ 2.5 Uc ~ y/nuc- 
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AE 0.4 



0.2- 

0.5 1 1.5 2 2.5 3 El. 

Figure 2: Shift function for n = 4 (solid) and its analytic approximation (^8[) (dashed). 
5.2 Quenching factors and infrared sensitivity 

The estimate (^) shows that in the presently available RHIC range p_L < 6 GeV a reliable 
quantitative prediction of quenching can be hardly made. It is the soft "singularity" of the LPM 
spectrum ujdl{uj) / duj oc 1/v^ that causes instability of the perturbative QCD description. 

At the same time, this very same instability makes the study of quenching in the p_L ~ 
20 GeV range, which is accessible at RHIC, the more interesting and valuable. In this region the 
characteristic gluon energies are comfortably large, tui ~ 2 GeV, while the quenching factor at 
the same time is still quite sensitive to much smaller energies deep into the infrared domain. 

To illustrate this point in Fig. ^ we show the expected quenching factors for the hot medium 
(g = 0.2 GeV^, sizes L = 2 and L = 5 fm) as a function of an infrared gluon energy cutoff. A 
sharp cutoff is not realistic and can be used only as a means of quantifying sensitivity of the 
answer to the low momentum region. In reality, what matters for quenching is the spectral 
properties of the (hot) medium which determine how the gluons with energies of the order of 
c^BH are being produced. 

The label "ra floating" in Fig. ^ (and below) means that for this p± range we have used the 
realistic fit to the vacuum spectrum provided by the PHENIX collaboration 0, 

J ^ vacuum ('„ N 

= const -(1.71 + pAGeV])-''-'\ 

dpj_ 

and evaluated the effective exponent n{p±) according to ([1|). 

In Fig. ^ the quenching factors are shown for large transverse momenta, where we have set 
n = 4, the asymptotic value. A mismatch between the Q values at the common point p± = 
20 GeV in Figs. |^ and ^ is due to the much stronger bias effect in the former case: n ~ 12 ^ 4. 
Knowing the p± dependence of the vacuum production cross section is, obviously, a necessary 
ingredient of a reliable quenching prediction. 
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Q n floating, L=2 fm Q n floating, L=5 fm 




Figure 3: "Infrared" dependence of the quenching factor for hot medium. The curves (from 
bottom to top) correspond to the gluon energy cuts 0, 100, 300 and 500 MeV. 

Comparing the two figures we remark that at large transverse momenta (and smaller n) the 
sensitivity to the infrared physics gets naturally reduced. 



Q n=4, L=2 fm Q n=4, L=5 fm 




Figure 4: Quenching factors for hot medium. (The curves are the same as in Fig. ^.) 

In Fig. ^ for the sake of comparison the magnitude of quenching for a final state matter 
with the transport coefficient equal that of the cold nuclear matter, q = gcow = 0.01 GeV'^, is 
displayed. 
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Q n floating, L=5 fm 



n=4, L=5 fm 





Figure 5: Quenching factors for cold medium. (The curves are the same as in Fig. ^) 

6 Conclusions 

In this paper we developed the perturbative QCD description of the suppression (quenching) 
of inclusive production of hadrons with large transverse momenta p± in a medium. We did 
not address the Cronin effect that enhances particle yield at moderately large p±. A simple 
analysis shows that parton rescattering effects both in the initial and in the final state (transverse 

momentum broadening) should die out as 0{qL/p\), to be compared with the O [l^/q/pl^ 
behaviour characteristic for quenching. 

We found that quenching is related with the multiplicity of medium induced (primary) gluons 
with energies larger than the characteristic energy uji as 



ln(5(p±) = / dze "^N^tuiz), uji{p±) 

'0 



dp. 



In 



dp] 



n > 1 . (49) 

n 



One of the results of this study is the observation that the "shift" parameter S in the commonly 
used parametrization (E^), 



da 



medium 



dpi 



d(T^^~(pX + 5) 

dpi 



equals neither the mean medium induced energy loss, S = AE oc L^, nor S = const ■ L. In fact, 
all over a broad kinematical region where quenching is sizeable, —lnQ[p±) = 0{1), the shift 
increases with p± as 



S{pa 



2asCF 



IT 



and equals twice a typical energy the quark looses to induced radiation. 
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We also observed that the particle (jet) quenching phenomenon is sensitive to the character of 
the distribution in the energy loss D{e) (in particular, to the energy region e ~ a'^oOc <^ oOc where 
D[e) has a maximum, see (^)), rather than to its high-energy tail e ~ tUc which determines the 
mean energy loss AE. The latter is dominated by a rare {O (««)) single gluon emission with 
maximal available energy, e ~ ~ cUc- 

Though formally a coUinear- and infrared-safe quantity, we found the quenching factor Q{p±) 
to be highly sensitive to the region of small gluon momenta, especially for p± < 20 GeV. The two 
effects that contribute to this are the bias (n ^ 1) and the LPM suppression which makes the 
gluon energy spectrum singular ai u 0. On the one hand, the region where the characteristic 
gluon energy Ui becomes comparable with cjbh ~ 300 MeV (which regulates the transition from 
the Landau-Pomeranchuk-Migdal to the Bethe-Heitler regime) the pure perturbative treatment 
is hardly applicable. On the other hand, the experimental (and possibly theoretical) studies 
of the momentum region corresponding to cui ~ (1 ^ 2) GeV which corresponds to p± = 10 ^ 
20 GeV (for n ~ 10) where the Cronin effect is expected to have disappeared, will provide an 
important information about the spectral properties of the final state medium produced in heavy 
ion collisions. 

The final remark concerns our treatment of a medium as uniform and static, which is obviously 
not a realistic approximation. For "hot" medium characterized by a large transport coefficient 
q ~ 0.2 GeV^ ^ 5 GeV fm~^ we estimate the characteristic lifetime of gluons with energies 
ui p±/n that determine quenching as 

^1 ~ \ ^— ~ I \ / '^—t: — 1 ffml . 

V"-^ V V 5-n GeVy ^ ^ 

In the most interesting and practically important situation this lifetime is small compared with 
the size of the medium, ti <^ L. In these circumstances it looks plausible that one will be 
able to calculate the quenching factor by employing the basic equations (pTsP with the dynamical 
position- dependent induced gluon distribution dl{u, z)/du!dz given in p8| . 
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